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Spin transport properties of the quantum one-dimensional non-linear sigma model are studied 
based upon the Bethe ansatz exact solution for the 0(3) sigma model and the 1/A'^-expansion 
approach for the 0{N) sigma model. It is shown that the spin transport at finite temperatures is 
ballistic and there is no spin difi^usion in this system. A possible mechanism for the spin diffusion 
observed in the NMR experiment of the Haldane gap system AgVP2S6 by Takigawa et al. [Phys. 
Rev. Lett. 76, 2173 (1996)] is also discussed. The inclusion of the spin-phonon interaction accounts 
for the experimental results qualitatively. 



I. INTRODUCTION 



D ' The spin diffusion in the Heisenberg spin chain is a long-standing issue in the studies of magnetisni.0 For the XXZ 
5^ 1 Heisenberg chajn-with spin s = 1/2, there is a strong piece of evidence for the absence of the spin diffusion in the 

■ massless state.B~Q In the case that the spin s is an integer, there exists an excitation gap, as predicted by Haldane many 
years ago using the mapping to the one-dimensional(lD) quantum 0(3) non-linear sigma model.Q Thus it is expected 

■ that the spin transport properties at low temperatures are quite different from those of the s = 1/2 Heisenberg chain. 
Recently, it was observed in the NMR experiment for the Haldane gap systerrL. AgVP2S6, by Takigawa et al. that 

I . a diffusive behavior appears in long wavelength spin-spin correlation functions. Q In order to explain the experiment, 
^ ' Sachdev and Damle showed that a semiclassical analysis of the quantum ID 0(3) non-linear sigma model gives a 
Q , diffusive behavior in the spin correlation function at low temperaturcs.EI However pthc ID quantum 0(3) sigma model 
^ ' is an integrable model, and exactly solvable in terms of the Bethe ansatz method.Ei The spin current is conserved, and 
thus it is expected that the spin diffusion is absent in this model. Although we can not calculate correlation functions 
directly by using the Bethe ansatz method, the spin transport properties are also probed by the Drude weight, namely, 
the zero frequency part of the spin conductivity. The calculation of the Drude weight for exactly solvable models was 
recently extended to the finite temperature case by Fujimoto and KawakamiJid If the Drude weight is non-zero at 
finite temperatures, the system shows a ballistic transport rather than a diffusive one. One of the purposes of this 
, paper is to discuss the exact Drude weight of the ID quantum 0(3) sigma model at finite temperatures based upon 
' the Bethe ansatz method. Our finding is that the Drude weight is non-vanishing at finite temperatures, and that 
\ the spin transport is ballistic. We will also confirm this result by utilizing the 1/A^-expansion analysis for the 0{N) 
non-linear sigma model. For large TV, we can calculate the spin correlation function perturbatively. We will show that 
the Ward-Takahashi identity for the current vertex operators leads to a ballistic form of the spin correlation function, 
"j^ ■ and that the spin diffusion does not exist. Thus in order to explain the spin diffusion observed in the Haldane gap 
I system AgVP2S6, we need to introduce a coupling with some other degrees of freedom, such as interchain couplings, 
i' , or spin-phonon interaction, etc. Moreover a remarkable feature of this experiment is the existence of a temperature- 
'■^ ■ dependent cutoff frequency, which can not be explained by the pure non-linear sigma model. The second purpose of 
^ , this paper is to investigate a possible mechanism of the spin diffusion in this system. The interchaig coupling J' of 
O ■ this material is so small, J' / J < 10^^, that it is unlikely to give rise to the observed spin diffusion.!] We will try to 
explain the experimental results by introducing the coupling with phonons. We will show that the inclusion of the 
spin-phonon interaction gives, at least qualitatively, a fine explanation for the spin diffusion and the presence of the 
temperature-dependent cutoff frequency observed in AgVP2S6. 

The organization of this paper is as follows. In Sec. II, the exact Drude weight at finite temperatures is derived based 
' upon the Bethe ansatz solution. In Sec. HI, the non-diffusive behavior of the spin correlation function is argued by 
using the 1/A^-expansion of the 0{N) sigma model. In Sec. IV, we discuss a possible mechanism for the spin diffusion 
observed in the NMR experiment taking into account a coupling with phonons. Summary and discussions are given 
in Sec.V. 
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II. THE BETHE ANSATZ EXACT SOLUTION FOR THE 0(3) SIGMA MODEL AND THE FINITE 

TEMPERATURE DRUDE WEIGHT 



In this section, we will derive the exact Drude weight of the quantum ID 0(3) non-linear sigma model at finite 
temperatures using the Bethe ansatz exact solution. The action of the 0(3) non-linear sigma model is given by, 

^ " ^ / d^dTt;(9^n)2, (1) 

where n — {rix^ny^nz) is a-unit vector, (n)^ = 1. The elementary excitations of this model are massive triplets 
interacting with each other .Ell The Bethe ansatz exact solution for this model was obtained by Wiegmann a decade 
ago.EI In order to obtain the Drude weight for the spin conductivit}?. we need to use the twisted boundary condition 
introducing the flux $ which couples with spin degrees of freedom£3'Ej The Bethe ansatz equations under the twisted 
boundary conditions are given by,0 



P=l 7 = 1 ' 

where p and q are defined by the eigen values of the Noether current operator I ~ nx dtfi/v; — p{p+ 1), Iz ^ q, 
and m cx eKp(—2Tr/g)/g. The total energy in the presence of a magnetic field h is given by 

E = my^cosh^a - h{p- q). (4) 

a 

The thermodynamic properties are obtained by using the string hypothesis for the scdjutions of the rapidities. As 
was shown by Wiegmann, 9a is real, and A^'*^ = '^S + ^T^in + 1 — 2fc)/2, fc = 1, 2, nn The Drude weight at finite 
temperatures is derived from the finite size corrections to the thermodynamic Betlie, ansatz equations. In order to 
calculate the finite size corrections, we expand the real part of the rapidities in l/iJlJ'tj 

Bg=eOg + J- + ^+ 0{1/L^), (5) 

AS = ASa + ^ + |^+0(l/L3). (6) 

Substituting these rapidities into eqs.(H) and (^), and taking a continuum limit, we obtain the coupled integral 
equations which determine the spectrum of the system. The zeroth order conttibutions in 1/L give the conventional 
thermodynamic Bethe ansatz equations which were first obtained by Tsvelik,ll3 

piB) + /(.) . - cosh. - J ^ d.'^^-^^p(.') 

nOO 

+ dAa„(0-A)a„(A), (7) 

a'^iX) + J2 ^nm * <Jrn{X) = / d0a„ (0 - A)p(0), (8) 

J — oo 

7^ 7ln(l + e-"^) 

_^ {B-B')^ + n^ 

(A) 

-TJ2 dAa„(0- A)ln(l-He-^), (9) 



m—l 



n=l 



ln(l + e^) = — - dBa„{B - A)ln(l -|- e~^) 



T 

^A„™*ln(l + e-^). (10) 
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Here p{9) and cr„(A) are the distribution functions for the rapidities O^a and Aq^, respectively, and p^{9) and (J^{\) 
are those for holes. n{9) = rin(p''/p), e„(A) = rin((T^/cr„), and 

2{n^l) 2(n+l) 

''"W 4^2 +^2(„_ 1)2 +4^2+^2(^+1)2' 



d f°° dx' X- x' 

Anm * 4>{x) ^ 5nm4l{x) + — J — 9„m( 



with 



^ 2 y V 2n - 2 / V 2n 

and 0(a;) = 2 tan^^ x. 

The finite size corrections to the above equations give the coupled equations which determine /i 2 and gi^2m 



r\ / \ I \ In — ml / \\n — ml + 2 J \n + m — 2) \n + m/ 

Qnm{x) = { H X 1^ X \ f X\ ^^^^ 



{p{e) + p\e))h{e) T d^^'Trr^iT^^— 2p(^')/i(^') 



2^ {d-e'Y + ^^' 

+ dAa„((?- AK(A)5i„(A), (12) 

n$ 

'^,'j(A)ffin(A) + ^ A„,„ * (t™(A).9i„(A) = / d0a„(0 - A)p(0)/i(0) , (13) 

™ — 1 J —00 ^ 



(P(^) + p''(0))/2(^) = - y ^ d^' (^^^7)7^P(^')/2(^') 

+ ^dAa„(0-A)a„(A)52„(A) + --[(p(0)+p"(0))/f(^?)] 



— 0' °° r°° 



da„(6' - A) 2 

cr"(A)gi„(A) 



dA 

+ (surface terms), (14) 
a,':(A)g2„(A) + ^ A„,„ * a™(A)52™(A) = / d0a„(0 - A)p(0)/2(0) 

rn=l 

+ ^^[('^(A)+a"(A))ffL]-E^*--5L+ r d0^^^i^p(^)/?(0) 

m J 00 

+ (surface terms). (15) 

The Drude weight at finite-|temperatures is given by the second derivative of the energy spectrum with respect to 
the Aharonov-Bohm flux 



2 \ d$2 



(16) 



$=0 



Here (• • •) is a thermal average. Note that the Drude weight D defined above is different from the second derivative 
of the free energy which is related to a Meissner fractionO i_. 
Evaluating the finite size corrections to the total energy up to the second order in l/ijl3 

E El E2 

we find that the dependence on $ appears in the second order term £'2- Then from eqs.(^) and (|7|)-(p^), we obtain 
the formula for the Drude weight, 
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D 



2 d$2 

^ /'OO 

2 



dK(6l) 



dA K(A)+a^J(A)) 



dffi„ 
d$ 



de„(A) 



This expression is similar to that for the Hubbard model 



dAVl + e^"W/^y a„(A) + (T4KA) dA 
.0 



(18) 



A. Low temperature expansion 

We now derive the temperature dependence of D in the low temperature limit. We consider the case of the non- 
vanishing magnetic field h ^ 0. In this case, from cq.([l^) we find e„(A) > 0. Thus, at low temperatures, eq.(||) is 
reduced to 



K{9)=mcoshe-h + j ^Ae'j^—^^^^K{9'), 



(19) 



where k{±B) = 0. In the case of /i < m, all excitations are massive, and the solution for the above equation is given 
by 

K{e) ^ mcoshe ~ h, (20) 

and B = 0. The low temperature limit of eq.([lO|) gives £ri(A) = nh and de„(A)/dA ^ exp(— (m — h)/T). From eqs.(|^) 
and (^), we have 

p[e) + p^{e) = ^cos\,e, (21) 

a„(A)+a:j(A)=e-^. (22) 
Moreover using eqs. (|l|) and (|l3|), we can show that at sufficiently low temperatures 

(23) 

K(A) + ^,':(A))%i^-^. (24) 
Substituting these expressions into cq.(|l8|), we obtain the Drude weight in the low temperature limit, 

D^Vfer"^. (25) 



Thus the Drude weight is non-zero at ffiiite temperatures. This result implies that the spin transport of the ID 
quantum 0(3) non-linear sigma model in the low temperature region is ballistic, and that the spin diffusion does not 
exist. 



B. High temperature expansion 

Here we consider the Drude weight in the high temperature limit. Although in the sufficiently high temperature 
region, the non-linear sigma model is not related to the s = 1 Heisenberg spin chain, it is interesting to check that 
the spin diffusion does not occur in the non- linear sigma model even in the high temperature limit. We mvoke the 
standard method for the high temperature expansion of the Bethe ansatz solution developed by Takahashi.Ej For this 
purpose, it is convenient to rewrite eqs. 

(0)> ®) and (^ into the following form. 
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p+p^ 


= — cosh 9 






= 52n3 * p 4 


-S*(CT^_l+Cr^+l), n 


CTi + a'l 


= S * CT2, 






= m cosh 9 - 


~ Ts * ln(l + e^J^), 


en(A) 


^Ts* ln(l 


+ e )(-!- + e ^ ) 



-d„,2rs * hi(l + e T^), ra=l,2,. 
Um — = h, 

n — >oo 77, 



(26) 

(27) 
(28) 
(29) 

(30) 
(31) 



where Eq = — oo, and 



f{x) = / da;' 



4cosh[7r(a; - x')/2] 



Since the non-hnear sigma model is a relativistic field theory model, it is appropriate to introduce the ultra-violet (UV) 
cutoff A which may be of order J, the exchange energy of the original s = 1 Heisenberg spin chain. In the case of 
T A, eqs.(E9) and (5|) are recasted into, 



lnC = -^ln(l + ry2), 
ln7?2 = i ln(l + 7?i)(l + ln(l + C'), 

\iirjn = + 7?n-i)(l + rin+i), n > 3, 

lim — — 0, 

n — >oo 71 



(32) 
(33) 
(34) 
(35) 



where C — / p, rjn — a^^/an- Thus in this limit, C and rjn are constants, and the solutions of the above equations are 
given by 



Vn = n" — 1 n > 2, 



(36) 
(37) 

Substituting these solutions into eqs.(|2^), (|27|), and (^), we can solve these coupled equations for p and cr„ by using 
the standard method. Then we find. 



Gttti 



2tt 



1 



1 



A)2 + (77_l)2 



A)2 + (77+1)2 



cosh 6', 



7i > 3, 



o'2(A) = 77;— s * cosh^ + 2s * (JsiX), 

cri(A) = S* (72 (A), 

p(0) = — -2s*a2(A). 



(38) 

(39) 
(40) 
(41) 



Here we have introduced the UV cutoff 9c = cosh "'^(A/77i) in order to suppress the divergence of the 0-integral in 
eq.(pa). We can also obtain d/i/d$ and d5i„/d<I> from eqs.(n2h and (O) in a similar manner. The results are 



d$ 



p 



5 

19^' 



0-1 ■ 
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d$ 



1 



1 

^ ~36^' 
77 > 3. 



0-2- 



d5 



12 



d$ 



1 

IStt' 



" d$ 37r77 ( 772 - 1) ' 

Furthermore in the high temperature limit, we have the following relations. 



(42) 
(43) 
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2tt dA 1 + e~ 
Then substituting these equations into eq.dlS), we have the Drude weight, 



1 Af~ 1 

p{e)^-sgnO--—^, (44) 
Ztt au 1 + ct 

/(^?) = ^sgn0^-^, (45) 
Ztt do 1 + eT 

a„(A) = ^sgnA^^-^, (46) 
a::(A) = :lsgnA^-^. (47) 



D-^, (48) 



with 



Note that D in the high temperature hmit has the same temperature dependence as that of the ID Hubbard model. 
Since even in the high temperature hmit, the Drude weight is non-zero, it is expected to be non- vanishing in the whole 
temperature region. Therefore the spin diffusion does not exist in the ID quantum 0(3) non-linear sigma model. This 
is the main conclusion of this section. 



III. 1/iV-EXPANSION APPROACH TO THE SPIN DYNAMICS OF THE 0{N) SIGMA MODEL 

In the previous section, we have calculated the Drude weight at finite temperatures which is a probe for transport 
properties. However in order to discuss the spin diffusion, it is desirable to study spin-spin correlation functions 
directly. This issue is not tractable in terms of the Bethe ansatz metJbod. Moreover, there is a criticism by Buragohain 
and Sachdev against the approach exploited in the previous section.ElJ They claimed that the calculation of the Drude 
weight from the finite size spectrum may not give the correct thermodynamic limit at finite temperatures. Thus, in 
this section, we will disc iiss . spin- spin correlation functions with the use of the 1/N expansion method for the quantum 
ID 0{N) sigma model,ElfE3 and check the consistency with the results obtained by thejJBethe ansatz method in Sec.H. 
The 1/A^-expansion is utilized by several authors in the study of Haldane gap systems. EJcJ Although the 0{N) sigma 
model with large N is not related to the Haldane gap systems, the qualitative properties of the excitation are similar 
to those of the 0(3) sigma model. It is expected that the results obtained for the 0{N) sigma model with large N 
may give a qualitative understanding of the 0(3) sigma model. 

We are mainly concerned with the spin lattice relaxation rate 1/Ti, which is given by 

-i-oclimVA^(,)I^^^, (50) 

111 uj^O^ — ' UJ 

q 

where x(g,a;) is the retarded spin correlation function, {S^{q,uj)S^{—q,—uj))^, and A{q) is a structure factor. It is 
known that in Haldana^ap systems, the main contribution to 1/Ti comes from the uniform part of spin-spin correlation 
functions for q ~ Ool£| It is given by the correlation function of the uniform magnetization vector I = n x dtn/v. 
That is, 

dx / dt([Z"(a;,i),r(0,0)])e'('^*+?^). (51) 

For sufficiently large N , the excitation in the 0{N) sijgina model is massive iV-plet bosons interacting weakly 
with each other via a coupling constant of order 0{\/N)]l:A Thus we can calculate the correlation function by the 
perturbative expansion in terms of The Green's function of the free massive A^-plet bosons n is given by, 

G^{p,LOn)^ r dT{Trn'^{q,r)n'^{-q,0)y^"^ ^ r ^2 f2 ' (52) 

Jo ['^^nV - 
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where oJn = 27rnT, and Ep = y'v'^p^ + rn?. One of the important effects of the interaction is the quasiparticle damping 
which is given by the imaginary part of the self-energy. The lowest order self-energy correction of which thej-daagram 
is shown in Fig.l, is easily evaluated as ImS^ ~ c{u>/T) exp(— 2rn/T) (c, a constant) for small iv and p ~ O.E3 If this 
damping effect is sufficiently large, it leads to the diffusive form of the correlator for the n field, 



1 



{uj + i(5)2 — i;2p2 _ jj)2 _ ic[uj/T) exp(— 2m/r) 



(53) 



Naively thinking, it seems that this correlation function gives rise to a diffusive behavior of 1/Ti. However this 
expectation is not correct. I is the conserved current of the non- linear sigma model which satisfies the continuity 
equation dtl + dxj = 0, where j = vn x dxti. Therefore a Ward-Tkahashi identity holds for these current operators. 
The Ward-Takahashi identity ensures the cancellation of the self-energy corrections with vertex corrections of the 
correlation functions, which leads to a ballistic behavior of spin-spin correlation functions. In order to show this, we 
introduce the three point vertex functions defined as. 



A^ip^,q^)G^ip^)G^q^^p^) = r{q^)n^{p^.Wiq^~p^)), 
A'}{p^,q^)G''{p^)G''{q,-p,) = {r{q^)nf\p^)n'^{q^^p,)), 

where = {jp,i£n) etc. In the momentum-energy space, the current operators are expressed as 



(54) 
(55) 



(56) 
(57) 



Using the commutation relations [P(x),nP(y)] — ie°'^"'n'^ {x)S{x — y), we can derive the Ward-Takahashi identity,Ea 

iw„A[^(p^,g^) +(jA^(p^,g^) = G'^q^, ^ p^,) - G-\p^). (58) 

Here we omit the spin index a of G'^{p^)^ since it does not depend on a by virtue of the 0{N) symmetry. We will take 
an analytic continuation iuJn — + -I- i^ at the last stage. Since we are concerned with the contribution from small q 
and u), s'n and p' in Aq i{p^, q^) can be replaced with e and p. As a result, we have A"i [p^, q^) « (u^p/ie„)AQ (p^, q^). 
Using this relation and the Ward-Takahashi identity (58), we can obtain the spin-spin correlation function. 



X{q ^ 0, iujm) = ^ Aq (p^, q^)G{p^)G{q^, ~ p^)i£^ 

Vu 



E 



lEr. 



IWm -I- 



.[G{p^)^G{q^-p^)]. 



Up to the lowest order in we expand G — Gq + Gq^Go in eq. (p9[) 
the limit of q — > and lu ^ 0, the above expression reduces to. 



Then, putting iLu„ 



(59) 



i(5, we find that in 



coth ^'^^ 



2T 



1 



16 



16 



-hregular terms. 



(60) 



Thus it has a non-diffusive pole at w = zL(v'^p/ Ep)q which can not be removed by finite temperature effects. This 
pole gives a logarithmic singularity ~ logw of l/Ti which signifies the ballistic jScin transport. The result is basically 
equivalent to that obtained for the free boson model by Jolicoeur and Golinclli.LJ This ballistic behavior is due to the 
cancellation between the damping effect and the vertex corrections. Although we have demonstrated the cancellation 
only up to the lowest order in 1/N, we expect that the cancellation holds in all orders. In this model the forward 
scattering is the dominant low-energy process. Thus the quasiparticle damping does not give rise a diffusive behavior, 
if one takes into account the vertex corrections which preserve the current conservation law. The above result is 
consistent with the Bethe ansatz solution obtained in Sec. II. Therefore we can conclude that the spin diffusion is 
absent in the ID quantum non-linear sigma model itself. 
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IV. COUPLING WITH PHONONS 



In this section, we will discuss about the origin of spin diffusion observed experimentally in the Haldane gap 
system, AgVP2S6. The important feature observed in this experiroent is that there exists the cutoff frequency oJc 
which depends on temperatures like ~ exp(r/r) with T a constant.Q At w > Wc, the spin diffusion is observed. One 
possible mechanism of the spin diffusion may be a coupling with some other degrees of freedom, such as phonons, as 
discussed by Narozhny for the s = 1/2 Heisenberg spin chain.Q Here we will investigate this possibility. In order to 
treat the interaction between spin and phonon, we invoke the 1/A^-expansion for the 0{N) sigma model. We add the 
following term to the Hamiltonian which represents the spin-phonon interaction, 

H' =g' f dx{b\x) + b{x)){d,n)\ (61) 



where b{x){b'' (x)) is an annihilation(creation) operator for phonons, and g' is a coupling constant. For the case of 
iV = 3, this term comes from the continuum limit of the spin-phonon interaction defined on lattices, 

H' ^J'Y,ibl + h)S,-S,+i. (62) 

i 

We expect that the analysis for the 0{N) sigma model gives a qualitative understanding for the effect of this term. 



For N oo, we can treat the spin-phonon interaction eq.(61) perturbatively. The diagram for the lowest order 
self-energy correction is shown in Fig. 2. We can easily evaluate this diagram using the Green's functions for massive 
magnons (|5^) and free phonons given by, 

2 2 

^b'-")- (— ^^7^' (63) 
where c is the velocity of phonons. Then the imaginary part of the retarded self-energy reads, 

ImS«(9 = 0,c.) ~ . ",„,.^rX' - C{T)a;, (64) 
Tsmh (§^) c** - 



for small ui. Here m* = raj — (v^ jc?). Then, if C(T) > Eq holds, the retarded Green's function of massive iV-plet 
n has a diffusive pole, 

G^\(l.^) = -. r, — 2T^' (65) 

where Ds — /C{T). In contrast to the case without phonons discussed in Sec. Ill, this damping effect does not 
cancel with vertex corrections of the spin-spin correlation function, eq. (pT|) , since the spin current operator is not 
conserved. Neglecting vertex corrections and using the diffusive propagator (^), we evaluate the RPA type diagrams 
for the spin correlation function (|5^). The result is given by, 

xiq,^^)^- 7^ , (66) 

for small q and uj. Here D'^, Uc, and Z have the same temperature-dependence as Ds- Thus we have a diffusive 
behavior of the spin-spin correlation function. The Fourier transformation of eq.(pq) with respect to ui leads to. 



f^^x(9,c.)e'"* = i?;e-[^=''^+-ll*l. (67) 

Therefore, tOc gives a cutoff frequency. For lu > lOc the diffusive behavior appears. The important point is that Wc 
depends on temperatures, uJc ~ Tsinh^(m*/2r), if we neglect the temperature dependence of the mass gap m. These 
results seem to give a qualitatively satisfying explanation for the spin diffusion and the temperature-dependent cutoff 
observed in AgVP2S6. Note that the presence of the temperature-dependent cutoff frequency is quite due to the mass 
gap of the spin excitations. The above argument is applicable only if the damping effect eq.(^J) is large enough to give 
rise to a diffusive pole. In this case, the perturbative approach may break down, and the corrections due to higher 
order processes are not negligible, though it is expected that the propagator has a diffusive form, eq. (|6^) , for small 
momentum and energy. Thus it may be difficult to fit our result to the experimental data precisely. However tuning 
the parameters, v, c, g' properly, we obtain the temperature dependence of ujc which is qualitatively consistent with 
the experimental results. 
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V. SUMMARY AND DISCUSSIONS 



In this paper, we have shown that the spin transport of the ID quantum non-hnear sigma model is balhstic even 
at finite temperatures, and that the spin diffusion does not exist. We have also argued a possible mechanism for the 
spin diffusion observed experimentally in the Haldane gap system AgVP2S6. It has been shown that the coupling 
with phonons accounts for the experimental results qualitatively. 

Our study is based upon the non-linear sigma model which is believed to be a low-energy effective theory for 
Haldane gap systems. We can not answer the question whether the spin diffusion intrinsically exists or not in the 
s = 1 Heisenberg spin chain model. According to de Alcantara Bonfim andJleiter, the ID classical Heisenberg spin 
model does not show the spin diffusion even in the high temperature limit .E3 The result for the classical spin system 
may be independent of whether the spin s is an integer or a half-integer, and generally quantum effects may suppress 
the thermal spin diffusion. Thus it is unlikely that the spin diffusion appears in the s = 1 quantum Heisenberg spin 
chain itself. We speculate that the spin diffusion observed in AgVP2S6 may be due to a coupling with some other 
degrees of freedom. However we need further investigations in order to reveal whether the spin diffusion is possible 
or not in the pure Heisenberg spin chains with an integer spin s. 
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FIG. 1. The diagram for the lowest order self-energy of massive A^-plct bosons. The solid line is the propagator for the 
massive AT-plet. The broken line represents the interaction between the A'^-plet bosons. 




FIG. 2. The diagram for the lowest order self-energy correction due to phonons. The broken line represents the phonon 
propagator. 
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